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Abstract 

We establish precise upper and lower bounds for the subelliptic heat kernel on nilpotent Lie groups G of 
H-type. Specifically, we show that there exist positive constants Ci, C2 and a polynomial correction function 
Qt on G such that 

CiQte"^ <Pt< C2(3te~^ 

where pt is the heat kernel, and d the Carnot-Carathcodory distance on G. We also obtain similar bounds 
on the norm of its subelliptic gradient |Vpt|. Along the way, we record explicit formulas for the distance 
function d and the subriemannian geodesies of H-type groups. 

On donne des estimations precises des bornes superieures et inferieures du noyau de la chaleur sous- 
elliptique sur les groupes de Lie nilpotents G de type H. Plus precisement, on montre qu'il existe des 
constantes positives Ci et C2, et une fonction polynomiale corrective Qt sur G telles que 

CiQte^^ <Pt< C2Qte~^, 

ou Pt est le noyau de la chaleur, et d est la distance de Carnot-Caratheodory sur G. On obtient aussi des 
estimations similaires pour la norme du gradient |Vpt|. En passant, on donne aussi des formules explicites 
pour la distance d et les geodesiques sous-riemannienes sur les groupes de type H. 
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1. Introduction 

Nilpotent Lie groups have long been of interest as a natural setting for the study of subelliptic operators; 
indeed, as shown in [24;], they model, at least locally, a general class of hypoelliptic operators on manifolds. 
Perhaps the simplest example is the classical Heisenberg group of dimension 3, followed by the higher- 
dimensional Heisenberg or Heisenberg- Weyl groups of dimension 2n+l having 1-dimensional centers. Beyond 
this, a natural generalization of the Heisenberg groups is given by the H-type (or Heisenberg-type) groups, 
which were introduced in [15]; these have a greater variety of possible dimensions while retaining some fairly 
strong algebraic structure. 

The main result of this paper is found in Corollary 14.31 in which we establish precise upper and lower 
pointwise estimates on the subelliptic heat kernel pt for an H-type group G, of the form 

GiQte-^ <pt< CaQte-^ (1.1) 
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for some positive constants Ci, C2 and an explicit function Qt, where d is the Carnot-Caratheodory distance 
on G. Additionahy, in Theorem l4.41 we obtain similar bounds for the subriemannian gradient of pi, namely 
that 

CiQ'e-T < iVpil < CaQ'e-T (1.2) 

for another explicit function Q' , where the inequality is valid at points sufficiently far from the identity of 
G. 

Estimates of the form p.ip for the classical Heisenberg group first appeared in in the context of 
a gradient estimate for the heat semigroup, as did an estimate equivalent to the upper bound in (|1.2p . A 
proof for Heisenberg groups in all dimensions followed in [20| . Our proof is similar in spirit to the latter, 
in that it relies on the analysis of an explicit formula for pt using steepest descent methods and elementary 
complex analysis. 

Less precise versions of the inequalities (jl.ip are known to hold in more general settings. Using Harnack 
inequalities one can show that for general nilpotcnt Lie groups, 

CiRi{t)e-^ <pt< C2{€)R2{t)e-T^ (1.3) 



for some constants c, Ci, C2 and functions i?2, where C2 depends on e > 0; see chapter IV of [27[. 
among others, improves the upper bound to 

Pt{9)<CR3{g,t)e-^, (1.4) 

with R a polynomial correction, using logarithmic Sobolev inequalities, whereas [iit improves the lower 
bound to 

Pt>G{e)R4(t)e'TT^. (1.5) 

Similar but slightly weaker estimates were shown for more general sum-of-squares operators satisfying 
Hormander's condition in [l7^ by means of Malliavin calculus, and in flij by more elementary methods 
involving homogeneity and the regular dependence of pt on t. 

In the specific case of the classical Heisenberg group, asymptotic results similar to (II. ip had been pre- 
viously obtained in [c?] and fn\, but without the necessary uniformity to translate them into pointwise 
estimates. A precise upper bound equivalent to that of (jl.ip was given in [Ij for Heisenberg groups of all 
dimensions. All three of these works, like j23| and the present article, were based on an explicit formula for pt 
and involved steepest descent type methods. In [s^, similar techniques were used to obtain a Li-Yau-Harnack 
inequality for the heat equation on Heisenberg groups. 

The proof we shall give here is largely self-contained, except for the formula (|4.2p for pt , which has been 
derived many times in the literature by many different techniques. We have also tried to err to the side of 
including relevant details. 

The author would like to extend his deep gratitude to his advisor, Bruce Driver, for his constant feedback, 
encouragement and support during the preparation of this paper. The author would also like to thank M. 
Salah Baouendi who provided some French linguistic advice, and the anonymous referee who suggested 
some useful references, especially Q. This research was supported in part by NSF Grants DMS-0504608 
and DMS-0804472, as weU as an NSF Graduate Research Fellowship. 



2. H-type groups 

H-type groups were first introduced in [isj. Chapter 18 of Q contains an extended development of their 
fundamental properties; we follow its definitions here, and refer the reader there for further details. 

Definition 2.1. Let g be a finite dimensional real Lie algebra with center 3 7^ 0. We say q is of H-type (or 
Heisenberg type) if Q is equipped with an inner product (•, •) such that: 



1- [3^,3^] 3 ; and 
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2. For each z G 3, define Jz ■ j ^ 3 by 

{Jzx,y) ^ {z,[x,y]) (2.1) 
where x,y E 3^. Then Jz is an orthogonal map whenever (z, z) — 1. 
An H-type group is a connected, simply connected Lie group whose Lie algebra is of H-type. 

Some authors use instead of item 2 the equivalent property that for a; G 3^ with ||a;|| = 1, the map 
a.dx '■ (kerada:)^ — > 3 is an isometric isomorphism. 

We record some algebraic properties of the maps Jz which will be useful later. We use \z\ :— ^ (z, z) to 
denote the norm associated to the inner product on g. The proofs are elementary and are omitted. 

Proposition 2.2. //g is a H-type Lie algebra, then the maps Jz '■ ^ 3^ defined in Definition \2.1\ enjoy 
the following properties: 

1. For each z, Jz is a well-defined linear map, and z i—^ Jz is also linear. 

2. — Jz • 

3. J^ — — ||z||2/. Thus for z ^ 0, Jz is invertible and J^^ — — llzH^'^J^. 
4- JzJw + JwJz= -^{z,w)l. 

5. {JzX,JniX)^{z,w)\\x\\'^. 

6. [x, Jzx] ~ WxW^z. 

Note that items [5] and [3] say that for z 0, Jz is an invertible skew-symmetric linear transformation of 
3-*-. Thus dim 3-*- must be even. We will write dim3'^ = 2n and dim 3 = m. 

Item m says that the subalgebra of End(3-'-) generated by the maps Jz is a Clifford algebra. In fact, 
it is a 2n-dimensional representation of C£o,m{^), the Clifford algebra generated by a real vector space of 
dimension m with a negative definite quadratic form (whose signature is (0, m)). So in order for an H-type 
algebra with dim3 = m, dim 3-'- — 2n to exist, it is necessary that C£o,m(lR) have such a representation. This 
condition is also sufficient: given such a representation, let V be the m-dimensional generating subspace of 
Cio,mO^), and let g = K^" ® V, with the maps Jz defined by the representation. Then the bracket on g can 
be recovered in terms of the Jz from 12.11 and g is an H-type Lie algebra. 

The Hurwitz-Radon-Eckmann theorem, as found in Q , gives necessary and sufficient conditions on n and 
m for such a representation to exist. The corresponding theorem for H-type algebras appears as Corollary 
1 of [1^, which we quote here. 

Theorem 2.3. For any nonnegative integer k, we can uniquely write k ~ a2^''+'? where a is odd and 
< q < 3; let p{k) := 8p 2^. (p is sometimes called the Hurwitz-Radon function.) There exists an H-type 
Lie algebra of dimension 2n -\- m with center of dimension m if and only if m < p{2n) . In particular, for 
every m G N there exists an H-type Lie algebra with center of dimension m. 

The special case m = 1 gives the so-called isotropic Heisenberg groups (also called the Heisenberg-Weyl 
groups) of real dimension 2ri + l; the very special case n = m = 1 is the classic Heisenberg group of dimension 
3. 

A Lie algebra g is said to be nilpotent of step fc if fc is the smallest integer such that all fc-fold brackets of 
elements of g vanish. A nilpotent Lie algebra is stratified if we can write g = gi® • • -ffig/c where [gi, g^-i] = gi 
and [gi,gfc] = 0. An H-type Lie algebra is obviously stratified nilpotent of step 2, with g^ = 3-*-, g2 = 3. 

We recall that given a nilpotent Lie algebra g, there exists a connected, simply connected Lie group 
G whose Lie algebra is g, and G is unique up to isomorphism. Indeed, we can, and will, take G to be g 
equipped with the group operation o given by the Baker-Campbell-Hausdorff formula, which for g nilpotent 
of step 2 reads 

xoy ■.= x + y+ -[x,y]. (2.2) 



In this case the exponential map g — > G is just the identity. It is obvious, then, that if are isomorphic 
as Lie algebras, then (fl,o), (g',o') as defined above are isomorphic as Lie groups. 

On the other hand, g can be identified as an inner product space with Euclidean space R2n+m ^ identifying 
^-^ with the first 2n coordinates and 3 with the last m. Therefore we can handle H-type groups concretely 
as follows. 

Proposition 2.4. If G is an H-type group, then there exist integers n,m and a bracket operation [•,•] 
on R2n+ni such that (R^""*"™, [•, •]) is an H-type Lie algebra whose center is R™ and G is isomorphic to 
(R2"+",o), where o is defined by i2.S\) . 

Henceforth we shall assume that any H-type group G is of this form. We shall use the notation 
g = {x, z) — (x^, . . . , x^", z^, . . . , 2™) to refer to points of G. The identity of G is (0, 0), and the inverse 
operation is given by {x,z)^^ — {—x,—z). Because of the identification of G with its Lie algebra, we will 
view [•, •] as a bracket on G. By a slight abuse of notation, we will also use [•, •] to refer to the restriction of 
[•, •] to R^" ® R^" C G ® G, which is a bihnear skew-symmetric mapping from R^" ® R^" to M™. The maps 
{Jz '■ z € R™} are identified with 2n x 2n skew-symmetric matrices which are orthogonal when \z\ = 1. 

We let {ei, . . . , e2n} denote the standard basis for R^", and {ui, . . . , Um} denote the standard basis for 



Note that the group operation on G does not preserve the inner product, and the vector space operations 
g i—f g -\- h, g eg are not group homomorphisms of G. However, the dilation 



is both a group and a Lie algebra automorphism for all a ^ 0. 

We can now identify g with the set of left-invariant vector fields on G, where Xi(0) — ^7, Zj{0) — ^j; 
then spanjXi, . . . , X2n} = 3^, spanjZi, . . . , Zm} — y We can compute 



R^ 



(pa{x,z) :— {ax,a^z) 



(2.3) 



(X,/)(x,z) = ^|t=o/((a;,^)o(ie„0)) 

= ^\t=of{x + tei,z + ^t[x, Ci]) 




So we have 




(2.4) 



(2.5) 



The (sub-) gradient on G is given in these coordinates by 




(2.6) 



Note in particular, if / is radial, so that f{x, z) ~ f{\x\ , \z\), this becomes 



V/(a;,z) = f\x\i\x\ , \z\)x + -f\z\i\x\ , |z|) \x\ J^x 



(2.7) 



where we use the notation u -A- to denote the unit vector in the u direction. We draw attention to the 

\u\ 

fact that X and JzX are orthogonal unit vectors in M^" for any nonzero a;, z. 
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3. Subriemannian geometry 



Our desired estimate for the heat kernel pt is in terms of the Carnot-Caratheodory distance d, which 
is best described in the language of subriemannian geometry. The goal of this section will be to obtain an 
explicit formula for d, and along the way we record formula for the geodesies of G. The computation is a 
straightforward application of Hamiltonian mechanics, but we have not seen it appear in the literature in 
the case of H-type groups. The corresponding computation for the Heisenberg groups (where the center has 
dimension m = 1) appeared in [2] as well as H: a computation for m < 7, which could be extended without 
great difficulty can be found in the preprint [4|]. 

Definition 3.1. A subriemmanian manifold is a smooth manifold Q together with a subbundle Ti. of TQ 
(the horizontal bundle or horizontal distribution, whose elements are horizontal vectors) and a metric (•, ■)q 
on each fiber Tig, depending smoothly on q G Q. Ti is bracket- generating at q if there is a local frame {Xi} 
for H near q such that spa,n{Xi{q), [Xi, Xj]{q), [Xi, [Xj,Xk]]{q), ■■■} = TqQ. 

An H-type group G can naturally be equipped as a subriemannian manifold, by letting Tig := {X{g) : 
X G 3^}, and using the inner product on g as the metric on Ti. In other words, Tig is spanned by 
{Xi(g), . . . , X2n{g)}, which give it an orthonormal basis. The bracket generating condition is obviously 
satisfied, since q — -^-^ ® b^,?^]- 

Definition 3.2. Let 7 : [0, 1] — > Q be an absolutely continuous path. We say 7 is horizontal if 7(t) €E ^^(t) 

for almost every t £ [0, 1]. In such a case we define the length of 7 as ^(7) := ^(7(t), 7(i))-Y(t) dt. The 
Carnot-Caratheodory distance d : Q x Q ^ [0, 00] is defined by 

d{qi,q2) mf{^(7) : 7(0) = 91,7(1) = 92,7 horizontal}. (3.1) 

Under the bracket generating condition, the Carnot-Caratheodory distance is well behaved. We refer the 
reader to Chapter 2 and Appendix D of [22] for proofs of the following two theorems. 

Tiieorem 3.3 (Chow). IfTi. is bracket generating and Q is connected, then any two points qi,q2 G Q are 
joined by a horizontal path whose length is finite. Thus d{qi,q2) < 00, and d is easily seen to be a distance 
function on Q. The topology induced by d is equal to the manifold topology for Q. 

Tiieorem 3.4. If Q is complete under the Carnot-Caratheodory distance d, then the infimum in the def- 
inition of d is achieved; that is, any two points 91 , 92 G Q o,fe joined by at least one shortest horizontal 
path. 

For an H-type group, we obtain the following explicit formula for the distance. Note that by its definition, 
d is left-invariant, i.e. d(g,h) — d{kg,kh), so it is sufficient to compute distance from the identity. By an 
abuse of notation, we write d{x, z) to mean rf((0, 0), {x, z)). 

Tiieorem 3.5. Define the function : K ^ M by 

^ 2^- sin 20 = ^ _ eot^^ = -^AOcoXe] (3.2) 
^ ' 1- cos 26* sin^e dO^ ^ ^ ' 

where the alternate form comes from the double-angle identities. Then 

fN^, z^O,x^O 
d{x,z)=<[x[, z = (3.3) 

[yi^, x = 



where 9 is the unique solution in [0,n) to v{9) — ^pj-. 
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We note that it is apparent from (|3.3p that we have the scaUng property 



d{(pa{x, z)) — ad{x, z) (3-4) 

with (fi as in (|2.3p . 

One way to compute the Carnot-Caratheodory distance is to find such a shortest path and compute its 
length. To find a shortest path, we use Hamiltonian mechanics, following Chapters 1 and 5 of [l^l- Roughly 
speaking, it can be shown that a length minimizing path also minimizes the energy ^ \\i{t)\\ dt, and as 
such should solve Hamilton's equations of motion. The argument uses the method of Lagrange multipliers, 
and requires that the endpoint map taking horizontal paths to their endpoints has a surjective differential. 
This always holds in the Riemannian setting, but is not generally true in subriemannian geometry; the 
Martinet distribution (see Chapter 3 of ,22]) is a counterexample in which some shortest paths do not satisfy 
Hamilton's equations. Additional assumptions on Ti. are needed. One which is sufficient (but certainly not 
necessary) is that the distribution be fat: 

Definition 3.6. Let 8 be the canonical f-form on the cotangent bundle T*Q, uj = dQ the canonical 
symplectic 2-form, and let := {pq £ T*Q : pq{'Hq) — 0} be the annihilator of H. (Note is a sub- 
bundle, and hence also a submanifold, of T*Q.) We say Ti, is fat if is symplectic away from the zero 
section. That is, if G TiP is not in the zero section, v G Tp^TiP , and u){v^w) = for all other w G Tp^Tf, 
then V — 0. 

Definition 3.7. If (Q,Ti., (•, •)) is a subriemannian manifold, the subriemannian Hamiltonian H : T*Q — > R 
is defined by 

H{pq) =Y,Pq{v.)^ (3.5) 

i 

where {vi} is an orthonormal basis for (Tiq, (•, ■)q). It is clear that this definition is independent of the chosen 
basis. Let the Hamiltonian vector field Xh on T*Q be the unique vector field satisfying dH + ll}{Xh, ■) = 
(as elements of T*T*Q). Xh is well defined because to is symplectic. Hamilton s equations of motion are 
the ODEs for the integral curves of Xh- 

The following theorem summarizes (a special case of) the argument of Chapters 1 and 5 of f23| . 

Theorem 3.8. IfTi. is fat, then any length minimizing path a : [0, f ] — > Q, when parametrized with constant 
speed, is also energy minimizing and is the projection onto Q of a path 7 : [0, 1] T*Q which satisfies 
Hamilton's equations of motion: 'y{t) — Xnijit))- 

We now verify explicitly that this theorem applies to H-type groups. We first adopt a coordinate system 
for the cotangent bundle T*G. 

Notation 3.9. Let (x, z, 77) : T*G M^n ^ M™ x M^" x M™ be the coordinate system on T*G such that 
x'ipg) = x'{g), 7J{pg) = zi{g), i,{pg) =p(^), rij{pg) =pi£j). That is, 

Pg = x{g),z{g), ^ ^,dx' + ^ 7]jdz^ 

\ i j 

In these coordinates, the canonical 2-form uj has the expression uj — d^i A dx'' + drjj A dz^ . 

Proposition 3.10. If G is an H-type group with horizontal distribution Ti. spanned by the vector fields Xi 
then TL is fat. 



■Hi C T;G by 



Proof. For an H-type group G, we have Pg E Ti.^ iff Pg{Xi{g)) = for all i. We can thus form a basis for 

= dz^ -"^dz^{Xi{g))dx' 

i 

" ^^■^ ^ X! ei) dx\ 



2 

i 
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Expressing pg in this basis as Pg = djW^ yields a system of coordinates (x, for Ti*^, where can 
be identified with the element {6^, . . . , 6*™) of R™. In terms of the coordinates (a;, z, ^, ?/) for T*G, we have 
r] = 9, C = -\Jex. 

So let 7 : (— e, e) TiP be a curve in which avoids the zero section. 7(0) is thus a generic element of 
TTi^ . We write ^{t) in coordinates as {x{t), z{t), 0{t)), where 9{t) 7^ 0. In terms of the coordinates {x, z, rf) 
on T*G, we have ri{t) — 9{t),£^{t) — ~^Jg(^t)x{t). Differentiating the latter gives 

Suppose that for all other such curves 7' with 7'(0) = 7(0), we have ^(7(0), 7'(0)) = 0. In terms of 
coordinates, 

= c.(7(o), 7'(o)) = ^(6(o)i'xo) - mi-m + J2ir,,mHo) - v-m'i^)) 
= m, x'm - (c'(o), i(o)) + im, i'm {v'm, m) 

= -^(-^^(0)2^(0) + Je(o)i(0),i'(0)) + ^(Je'(o)^'(O) + ^fl'(o)i'(0),i(0)) 

+ {m,z'm~{9'io),m) 
+ m,z'{o))-{9'io),m) 

For arbitrary u £ M™, take j'{t) = (a;(0),z(0) + tu,9{0)); then = ^(7(0), 7'(0)) = {9{0),u), so 
we must have ^(0) = 0. Next, for arbitrary v S R^", take 7'(i) = (a;(0) + tv, z{0),9{0)); then we have 
= ( J5i(o)U, i(0)). But 9(0) 7^ by assumption, so Jg^, is nonsingular and we must have i(0) = 0. 
Finally, take Y{t) = (a::(0), z(0), 6'(0) + tu); then (u, i(0)) = 0, so i(0) = 0. Thus we have shown that if 
Ci;(7(0), 7'(0)) = for all 7', we must have 7(0) ~ 0, which completes the proof. □ 

We now proceed to compute and solve Hamilton's equations of motion for an H-type group G. 
The subriemannian Hamiltonian on T*G is defined by (c.f. p.Sp ) 



^ 2n 



e t;g. 



(3.6) 



In terms of the above coordinates, we may compute 



Pg{Xi{g)) ^pg 




so that 



Recall that a path 7 : [0, T] T*Q satisfies Hamilton's equations iff j{t) = Xnijit)), i.e. dH^^j.-^ + 

^(tW,-) = 0. 

In an iJ-type group G, we write 7 in coordinates as "f{t) — {x{t), z{t),^{t),r]{t)) : [0,T] T*G, so that 
we have 



Thus Hamilton's equations of motion read 



.,_dH_ x__dH_ . . dH 



(3.7) 



To compute the derivatives, we note that ^V^ \Ax + = A* Ax + A*y. If we write Bxi] — JrjX, then 
(B,r,,y) = (?7, so B* = [x,-], and B^B, = \x\^ I. So for a path 7(i) = {x{t), z{t),m,v{t)) ■ [0,T] -> 

T*G, Hamilton's equations of motion read 



X — \7 = + -JrjX 



1 2 1 

= 7 l^^l '7+ ^[a:,^] 



7) = -V^B = 0. 



(3.8) 

(3.9) 

(3.10) 
(3.11) 



Theorem 3.11. {x{t), z{t)) is the projection of a solution to Hamilton's equations with x(0) ~ 0, z(0) = 
and x{l), z(l) given if and only if: 



1. If z{\) — 0, we have 

2. If z{l) ^ 0, we have 



v{t) = tx{l), z{t) = 0. 



(3.12) 



where, if x{l) ^ Q we have 



where 9 is a solution to 



and if x{l) — we have 



for some integer fc > 1. 



z{t) 



770 = 29 



2 hoi' 
z(l) 



(|77o|i- sin(|77o|i))'?o 



Wi)l 



v{9) 



4|^(1)| . 

wi)r 



770 = 27rfc 



z(l)_ 
l<l)l 



(3.13) 
(3.14) 

(3.15) 
(3.16) 

(3.17) 



Proo/. We solve (I3IHH3IIII) , assuming x(0) = 0, z(0) = 0. By (fXTT|) we have r]{t) = 77(0) = 770- H f/o = 0, we 
can see by inspection that the solution is 



77(t)=0, ^(0=^0, x{t)^t^o, z{t)^Q, 



(3.18) 



namely, a straight line from the origin, whose length is clearly This is (|3.12p . which we shall see is 

forced when z(l) — 0. 

Otherwise, assume ijo ^ 0. We may solve p.Sp for ^ to see that 



C = i - - Jr^oX. (3.19) 



Notice that substituting p.lPp into p.9p shows that 



(3.20) 



from which an easy computation verifies that (a;(t), z(t)) is indeed a horizontal path. 
Substituting p.lOp into the right side of (|3.10p shows that 



since J^gX = — |?7o| x. Thus 

e=ij,,,x + ^o (3.21) 

where ^0 = C(0)- If — 0, it is easily seen that we have the trivial solution x{t) — 0, z{t) — 0, ^(t) — 0, 
?7(t) = 770, so we assume now that 7^ 0. p.2ip may be substituted back into (|3.8p to get 

x^J^gX + ^o (3.22) 

so that 

l^ol 

Differentiation (or substitution) shows 

X = e"™^o- (3.24) 

Note that 

1^1' = A - %r = A(l - cos(ho| 0) ICol' . (3.25) 

It is easy to see from (I3.23P that x(t) lies in the plane spanned by ^0 a-nd Jrio£,Oi ^nd x(t) sweeps out a 
circle centered at j^^JrjoS,o and passing through the origin. In particular, the radius of the circle is l^ol/hol- 



X = (J,J"i(e*^"o - 7)^0 - -^J,o(e*^™ - 1)^0- (3.23) 



Now substituting ((3?23| and ([3?24| into ([3?20| . we have 

1 



ICol 



2|% 

|2 

-2 (1 - cos(|77o|t))r?o 



2|%| 



By integration. 



In particular. 



'' = ^^i\m\t-smi\7^o\t))vo. (3.26) 
2 hoi 

kl = :774(l^o|t-sin(holO)- (3.27) 
2|%| 



We note that inspection of (|3.27p shows that z(t) 7^ for t > 0. Thus the only solution with z(l) = is 
that of ((XT^ . 

To make more sense of this, let r = |^o| / Ivol be the radius of the arc swept out by x(t), and 4> = |?7o| t 
be the angle subtended by the arc. Then 



1 



which is the area of the region between an arc of radius r subtending an angle cj) and the chord which spans 
it. 

We must determine ^0:'7o in terms of a;(l),z(l). We have already ruled out the case z(l) = 0. If 
x(l) = 0, then (I3.25|) shows we must have I770I = 2fc7r for some integer fc > 1. p. 261 I3.27P then shows 
?7o — 2fc7rz(l)/ 1^(1)1, and |^o| = \/4fc7r |2:(1)|, as desired. In this case the direction of ^0 is not determined 
and ^0 niay be any vector with the given length. 

On the other hand, if x(l) 7^ 0, then |r/o| is not an integer multiple of 27r, so we may divide p.27p by 
(P?^ to obtain 

^ = - -AO) (3.28) 

\x(lf 4(l-cos|ryo|) 4 ^ ^ 

taking 9 — ^ |77o|, where v is as in p.2p . Then by (|3.25p we have 

l^ol' = \ 1^(1)1' 1 h = 1^(1)1' (3-29) 

2 1 — cos(|77o|) sm^ 

Note that once the magnitudes of ?7o; Co a-re known, their directions are determined: r/o = ^^(1) |f]o| / 
by (|3.26p . while can be recovered from (I3.23P : 

i^^-vl{Jn,{e'^^' -I))-'x{l). 

So ?7o,Co and hence x{t),z{t) are all determined by a choice of I770I satisfying (|3.28p . Writing 9 = I770I gives 
(13T3J13T4]). 

The "if" direction of the theorem requires verifying that the given formulas in fact satisfy Hamilton's 
equations, which is routine. □ 

To prove Theorem l3.5[ we must now decide which of the solutions given in Theorem l3.11l is the shortest, 
and compute its length. We collect, for future reference, some facts about the function v of p.2p . 

Lemma 3.12. There is a constant c > such that v' {9) > c for all 9 G [0, tt). 

Proof. By direct computation, ^'{9) = ^^'''"f^^/g'^"^^'* ■ By Taylor expansion of the numerator and denominator 
we have i^'(O) — 2/3 > 0. For all 9 e (0,7r) we have sin'^ 6* > 0, so it suffices to consider y{9) :— sin9 — 9cos9. 
Now y(0) = and y'{9) = 6* sin 61 > for 61 e (0,7r), so y{9) > for 61 e (0,7r). Thus v'{9) > for 6* G [0,7r), 
and continuity and the fact that lime|„ ^'(0) = +00 establishes the existence of the constant c. □ 

Corollary 3.13. 1/(9) > c9 for all 9 £ [0, tt), where c is the constant from Lemma \3.1S\ . 

Proof. Integrate the inequality in Lemma [3. 121 Note that z/(0) = 0. □ 

Proof of Theorem \3.5\ We compute the lengths of the paths given in Lemma 13.111 The z = case is 
obvious. Observe that for a horizontal path a{t) = {x{t),z{t))^ we have &{t) = Y^=l^^{^)^i{l{^))^ so that 
||(T(t)|| — \x{t)\. For paths solving Hamilton's equations, p.24p shows that \x{t)\ = j^ol: so ^(7) = j^ol- In 
the case x = 0, we have j^ol — \/ ^kir \z(l)\, where k may be any positive integer; clearly this is minimized 
by taking k — 1. 

Now we must handle the case x ^ 0, z ^ 0. In this case we have £{-/) = |Co| = ^j^, by p.29p . where 9 
solves p.l7p (recall 9 = ^ \vo\)- The function has j^(0) — 0, 1^(71) = +00, and by Lemma [3. 121 v is strictly 
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increasing on [0, tt). Thus among the solutions of (|3.17p there is exactly one in [0, tt). We show this is the 
solution that minimizes (^j^) and hence also minimizes ^(7). 

For brevity, let y = -pd. If y e [0, 7r/2] then y ~ v{9) for a unique € [0,oo). This is because 
v{9) > i'{tt/2) = tt/2 for 6 > 7r/2. Since 9 is increasing on [0, tt) it suffices to show this for 9 > tt. But for 
such 9 we have 

61 - sin 61 cos 61 9-^ ^ 1 1 tt 

= —n >-^>0-n>^-n> 

sm 9 sm 9 



2 2 2 



since sin 9 cos 6* < i for all 



Otherwise, suppose y > 7r/2. Let 



F{9) :- 



(9) 61 -sin 61 cos 6* 

which is smooth on (7r/2,cx)) after removing the removable singularities. We will show that if 7r/2 < 9i < 
TT < 92, then F{9i) < F{92). Thus if 9i is the unique solution to y = v{9) in (7r/2,7r) and > tt is another 
solution, we will have 

1. \ 2 

v{0i)F{9i) = yF{9i) < yF{92) - v{92)F{92) 



^ sm Wi / \ sm 02 

Toward this end, we compute 

26'(6l - sin 61 cos 6*) - 6*2(1 - cos^ 9 + sin^ 9) 



F'{9) 



20 cos 9 (9 COS. 9 ~ sini 



For 9 e (7r/2,7r) we have cos9 < 0, sin 6* > and thus F'{9) > 0. So F{0i) < F{tt) and it suffices to show 
F{-k) — tt < F{02)- We have F'{-k) = 2 > so this is true for 02 near tt, and F{+co) — +00 so it is also true 
for large 02- To complete the argument we show that it holds at critical points of F. Suppose F'{0c) = 
where 0c > tt; then either cos 6'c = or 0c cos 0c — sm0c = 0. If the former then F{0c) = 0c > tt- If the 
latter, then 0c — tan 0c, so 

fl2 n2 ra 
plQ\ ^ ^ 0_c ^ 0c > > TT 

^ 6'c — sin 6'c cos 6'c 6'c — tan 6'c cos^ 6'c 6'c( 1 — cos^ 6'c) ~ 
which completes the proof. □ 

Notation 3.14. If /, ft, : G ^ M, we write /(g) x ft(g) to mean there exist finite positive constants Ci, C2 
such that Cih{g) < f{g) < C2h(g) for all g G G, or some specified subset thereof. 

Corollary 3.15. d{x,z) x \x\ + \z\^^^ . Equivalently, d{x,zf' x + \z\. 

Proof. By continuity we can assume a; 7^ 0, z 7^ 0. If is the unique solution in [0,7r) to v{0) ~ we 
have d{x, z)^ = \x\^ (il^) > if I'^t 

it will be enough to show there exist Di,D2 with < Di < F{0) < D2 for all E [0,tt). F is obviously 
continuous and positive on (0, tt). We can simplify F as 

m 



sin^ I 
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from which it is obvious that hmgiTr ^(^) = tt > 0, and easy to compute that hmgjo ^(^) = 1 > 0, which is 
sufficient to estabhsh the coroUary. 

□ 

Results of this form apply to general stratified Lie groups. A standard argument, paraphrased from 
where many more details can be found, is as follows. Once it is known that d generates the Euclidean topology 

1/2 

on G, then d{x, z) is a continuous function which is positive except at (0, 0). d'{x, z) :— another 
such function, so the conclusion obviously holds on the unit sphere of d! . Now d' {(pa{x^ '^)) — ad' (x^ z), and 
inspection of (|3.3p shows that the same holds for d, so for general (x, z) it suffices to apply the previous 
statement with a — d'{x,z)^^. 

4. The sublaplacian and heat kernel estimates 

Definition 4.1. The sublaplacian L for G is the operator given by 

i = E^' (4.1) 

i 

where Xi are as given in (|2.4p . The heat kernel pt for G is the unique fundamental solution to the corre- 
sponding heat equation (L— ■§f)u — 0; that is, pt — e*^(5o, where is the Dirac delta distribution supported 
at 0. 

L is obviously left-invariant. L is not strictly elliptic at any point of G, but it is subelliptic everywhere. 

If we view the left-invariant vector fields {Xi} as elements of the Lie algebra g of G, they are an 
orthonormal basis for j-*-, which generates g: that is, span{Xi, [Xj, X^] : i,j,k — 1, . . . , 2n} = g. (It is easy 
to see that L does not actually depend on the choice of orthonormal basis {X^} for j-*-, but only on the 
inner product (•, •) on g.) We thus have span{Xj(g), [Xj,Xk]{g) ■ i,j,k = 1, . . . ,2n} = TgG for each g <E G 
(it is obvious for g = 0, and for other g it follows by left invariance). Thus the collection of vector fields 
{Xi} is bracket generating. By a famous theorem of Hormander ([lO|), L is hypoelliptic; that is, if Lu is 
on some open set, then so is u. Another case of Hormander's theorem applies to the operator L — ^ 
on G X (0,oo) = {{g,t)}] thus, since (L - ^)pt = is G°°, pt itself is G~ on G x (0,cx)). 

Our next step is to record an explicit formula for pt{x, z). Various derivations of this formula appear in 
the literature. For general step 2 nilpotent groups, derived such a formula probabilistically from a formula 
in [3] regarding the Levy area process. Another common approach, worked out in [Hi, involves expressing 
Pt as the Fourier transform of the Mehler kernel. [2^ has a similar computation, [ij] obtains the formula 
for H-type groups as the Radon transform of the heat kernel for the Heisenberg group. Other approaches 
have involved complex Hamiltonian mechanics ([H), magnetic field heat kernels ([lf|), and approximation 
of Brownian motion by random walks ([12]). In our notation, we find that 

p*(x,z) = (27r)-"(4^)-" / e»(A,.)-i|A|coth(t|A|)N^ f }^\ Y dX. (4.2) 

Jr". Vsinh(i|A|)y 

We can see directly by making the change of variables X = a^X' (among other means) that 

pt{x, z) = a^^^^+"^p^2t{ax, a^z) = a2(™+")p„2,(^„(a;, z)). (4.3) 

In particular, taking a — i^^/^, 

pt(x,z)) =t-'"->i(ii/2:c,tz) =i-™->i(^,V2(a:,z)). (4.4) 

Therefore an estimate on pi will immediately give an estimate on pt for all t, and we study pi from this 
point onward. 
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We immediately note that the integrand in (|4.2p has even real part and odd imaginary part, so that pi 
is indeed real. Moreover, being the Fourier transform of a radial function, pi is radial, i.e. pi{x, z) depends 
only on |x| and \z\. So we can apply (|2.7p and differentiate under the integral sign to get 

Vpi(x, z) - -i(27r)-™(4^)-" \x\ {qi{x, z)x + q2{x, z)Jix) (4.5) 

where 

gi(x,.) = -n%rT^- / e»<^-^)-^H--l^ll^l^f J-y^'cosh(|A|)dA (4.6) 
\x\ d\x\ 7r™ Vsinh(|A|)y 

g2(x,z) = %^= / e^(^-)-il^l-*'>I^IN^f-}^)"(-z)(A,z)dA (4.7) 
d\z\ 7m™ Vsinh(|A|)y 

As before, (|4.6p and (|4.7p do not really depend on z but only on |a;| , |z|. 

We now state the main theorem of this paper: the precise estimates on pt and its gradient. The proofs 
will occupy the remainder of the paper. 

Theorem 4.2. There exists do > such that 

1 + {\x\ d[x, z))" 2 

for d{x, z) > do. 
Corollary 4.3. 

Pt(x,z)^t " " ^ ^J-^li — e it^^^^^-i (4.9) 

l + (i|x|d(a;,z))«"^ 

/or (x, z) €z G, t > 0, with the implicit constants independent oft as well as (x, z). 

Proof. Theorem 14.21 establishes (|4.9p for t = 1 and d{x,z) > do. For d{x,z) < do the estimate follows from 
continuity and the fact that pt{x,z) > 0. Although the positivity of pt is not obvious from inspection of 
()4.2p . it is well known. A proof of this fact could be assembled from the fact that the semigroup e*^ is 
positive and hence pt > (see, for instance. Theorem 5.1 of 13]) together with a Harnack inequality such 
as Theorem III.2.1 of 27] (which is written about positive functions but easily extends to cover those which 
are nonnegative). 

Once ()4.9p holds for all {x,z) and t — 1, (|4.3p and p.4p show that it holds for all t, with the same 
constants. □ 

We also obtain precise upper and lower estimates on the gradient of the heat kernel. Again we work only 
on d{x, z) > do, and since Vpt vanishes for a; = 0, it is not as clear how to extend to all of G. However, the 
upper bound is sufficient to establish (14. lip , which is of interest itself. 

Theorem 4.4. There exists do > such that 

J/^ ^2n— m+l 

\Vpiix,z)\-\x\ ^ ' ' (4.10) 

' ^ ^' ' ^ l + {\x\d{x,z)r+'^ ^ ^ 

for d{x,z) > do. In particular, we can combine this with the lower bound of Theorem \4.^ to see that there 
exists C > such that 

lVpi(x, z)\ < C(l + d{x, z))pi{x, z). (4.11) 

The function q2 is of interest in its own right, because it gives the norm of the "vertical gradient" of 
Pi- I92I = {{ZiPi, . . . , ZmPi)\. The proof of Theorem 14.41 includes estimates on 52; we record here the upper 
bound. 

13 



Theorem 4.5. There exists do > and a constant C > such that 

\{Zm, . . . , Zr.p^){x, z)\ = Mx, z)\ < C )\> . (4.12) 

1 + (\x\ d{x, z)y^ 2 

whenever d(x, z) > da. In particular, for all {x, z) G G we have 

\{Zipi,...,ZmPi){x,z)\ < Cp{x,z). (4-13) 

Remark. Since our estimate is based on analysis of the formula l|4.2p . we will henceforth treat (|4.2p as the 
definition of a function pi on R^"+™. In particular, it makes sense for all n,m, whether or not an H-type 
group of the corresponding dimension actually exists (which can be ascertained via Theorem 12. 3p . The 
proofs of Theorems 14.21 and 14.41 do not depend on the values of n and to, so they likewise remain valid for 
all n, m. The estimates given are in terms of the distance function d, which likewise should be taken as a 
function defined by the formula (|3.3p . Indeed, the only place where we need pi to be a heat kernel is in the 
proof of Corollarv l4.3( where we use the positivity of pi which follows from the general theory. 

In particular, in Section [7] we shall make use of estimates on pi for values of n, to not necessarily 
corresponding to H-type groups. 

The proofs of these two theorems are broken into two cases, depending on the relative sizes of and 
\z\. Section E] deals with the case when \z\ < \x\ ; here we apply a steepest descent type argument to 
approximate the desired function by a Gaussian. Section E] handles the case \z\ 3> |a;|^ by a transformation 
to polar coordinates and a residue computation which only works for odd to. The result for to even can be 
deduced from that for to odd by a Hadamard descent approach, which is contained in Section [71 



5. Steepest descent 

We first handle the region where \z\ < Bi \x\^ for some constant Bi. li 9 = d{x, z) is as in Theorem 13.51 
this implies i'{9) < 4i?i; since v increases on [0,7r) we have < 9 < 9q in this region. Note also that by 
CoroUarv 13 . 1 51 we have d{x, z)"^ < £'2(1 + Bi) \x\ , as well as d{x, z)^ > |x| which is clear from p.3p . Thus 
for this region the bounds of Theorems 14. 2[ 14.41 and 14.51 are implied by the following: 

Theorem 5.1. For each constant Bi > there exists do > such that 

pi(x,z)x^e-^''(-^^)' 
\x\ 

h{x,z)\<-^e-i''^--^^'>\ i = l,2 
\x\ 

C 1 2 

_L.f.-\d(x,zf ^ niax{|gi(a;, z)\ , \q2{x, z)\} 
\x\ 

for all X, z with d{x, z) > do and \z\ < Bi \x\^ . 

Our approach here will be a steepest descent argument. Very informally, the motivation is as follows: 
given a function F(a:) = J^e'^ ^'■^^a{X)dX, move the contour of integration to a new contour T which passes 
through a critical point Ac of /, so that /(A) « /(Ac) + ■^f"{Xc){X — Xc)^. Then we have 

F{x) « e--'^(^=) ^ e--'^"(^=)(^-^=)'/2a(A) dA. 

For large x the integrand looks like a Gaussian concentrated near Ac, so F(x) x e~^ ^^^"'^ — "^'^"•^ Our 

proof essentially follows this line, in instead of R, but more care is required to establish the desired 
uniformity. 

14 



(5.1) 
(5.2) 
(5.3) 



Our first task is to extend the integrand to a meromorphic function on C™, so that we may justify 
moving the contour of integration. 

Let • denote the bilinear (not sesquilinear) dot product on C™, and for A £ C™ write := A • A; this 
defines an analytic function from C™ to C, and A^ = |A| iff A G K™. For w € C, let y/w denote the branch 
of the square root function satisfying Im ^/w > and ^/w > for w > (so the branch cut is the positive 
real axis). Thus if g : C — s- C is an analytic even function, A i— > g{y/)^) is analytic as well, and satisfies 
5(^/A2) = 5(|A|) for A e M". This holds in particular for the function ^is^, and thus the functions — 

and VA^coth VA^ are analytic away from points with = ikn, fc = 1, 2, . . . . 
Using this notation, we let 

VA2 \ 



, \ n+1 



Y sinh V A-^ 

ai(A) := cosh Va^ 

\ sinh V A^ / 



As mentioned previously, z may be any unit vector in R™ without affecting the computation. Therefore we 
shall treat it as fixed, while |z| is allowed to vary. 

Also, for A e C™, 61 e [0, 9a], ze S""-^ C M", we define 



/(A, 0, z) := -iiy{0)X -z + VX^ coth V A2 (5.4) 

so that ^ 

^/(A,6'(a;, z), ^)^-iX.z + ^ VA^ coth VA^ . 
4 \z\ 4 

We henceforth write 9 for 9(x, z). Thus we now have 

pi(x,z) = (4^)-"-" / e^^^(^''''^)ao(A)dA (5.5) 
g,(a;,z) = (47r)-"'-" / e-^/(^'''^^'a,(A) dA, i ^ 1,2 (5.6) 



Written thus, the integrands have obvious meromorphic extensions to A G C", analytic away from the set 
{VA2 =i/c7r, A: = 1,2, ...}. 

A simple calculation verifies that ^locothui = ii'{~iw), so we can compute the gradient of / with 
respect to A as 

Va/(A, 9, z) = -iv(9)z + w{-i\f}?)X (5.7) 
which vanishes when A = i9z. Thus i9z is the desired critical point. We observe that 

fii9z, 9, z) = 91^(9) + i9 coth{i9) = 9(y{9) + coi{9)) = — 5- (5.8) 

sin 9 

so by (p^j) . 

\x\^f{i9z,9,z)=d{x,zf. (5.9) 



Thus we define 



02 

iP{X, 9, z) := /(A, 9, z) - f{i9z, 9, z) = -w[9)X ■ z + ^fx? coth ^ 5-. (5.10) 

sin 9 
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We then have 



and analogous formulas for gi , (72 • Thus let 



(5.11) 



(5.12) 



It will now suffice to estimate hi. 

The first step in the steepest descent method is to move the "contour" of integration to pass through 
iOz. Some preliminary computations are in order. 



Lemma 5.2. For a,b E 



have 



|a| - |6| < Re V(a+W < |a| , < Im ^/(a + bi)^ < \b\ . 
Equality holds in the upper bounds if and only if a and b are parallel, i.e. a = rb for some r G 
Proof. First note that (a + bi)^ = \a\^ — \b\^ + 2ia ■ b. So by the Cauchy-Schwarz inequality, 

\{a + bif\^ = -\b\^f + (2a-bf 
<{\a\'-\b\y+A\a\'\b\' 
= (l«l' + |6|')^ 



(5.13) 



(5.14) 



so that I (a + 6i)^| < |a|^ + |6|^. Equality holds in the Cauchy-Schwartz inequality iff a and b are parallel. 
On the other hand, 



\{a + biY\ > Re{a + biy = \a\^ - \b\' 



(5.15) 



Now we can write 



{Rey'{a + bi)^y = i (^^i^^+W + Via + Wy 



- i^{a + bif + {a + biY + 2 V(a + bif 
].{\a\'-\bf + Ua + Uf\). 



The upper bound for 
otherwise we have by (|5.15p that 



Re ^{a + biY then follows from (|5.14p . The lower bound is trivial if \a\ < \b\, and 



(Rc ^{a + bi)^y > \af - \bf > {\a\ 



The lower bound for Im ^yJa + bi)^ holds by our definition of ^/^, and the upper bound is similar to the 
previous one. □ 

Lemma 5.3. For each Oq e [0,7r) there exists c{9q) > such that if a,b E M" with \a\ > c{9o), \b\ < 2tt, we 
have 



and 



for all e e [0, 6*0], z,xe C 



Reip{a + ib,e,z) > \a\ /2 

\ai{a + ib)\ < 1 
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(5.16) 



(5.17) 



Proof. Fix 6*0 G [0,7r). Note first that 



Re ipia + ib, 9, z) ^ v{e)b • z - Re f{i6z, 9, z) + Re yj{a + hiY coth yj{a + bi) 



(5.18) 



By continuity, i^{9)b- z — Ke f{i9z, 9, z) is bounded below by some constant independent of a for all 9 £ [0, 9q] 
\b\ < 2tt. Thus it suffices to show that for sufficiently large |a|, 



Re 



V(a + 5i)2 coth 1/(0 + biy 



Now for a e M, /3 e [-27r, 27r] we have 



Re((a + iP) coth(a + = 



a sinh a cosh a + (3 sin /3 cos /3 



cosh a — cos^ P 



> a coth a 



> a coth ( 
>4H 



/3 



cosh a 



2tt 



cosh a 



for sufficiently large \a\. (Recall that limQ_^±oo cotha = ±1-) Thus, since 



Re y/{a + bi)^ > \a\ - \b\ > \a\ - 2-k 



(5.19) 



and 



Im ^(a + bi)'' 



< 2n, 



it is clear that (|5.19p holds for sufficiently large \a\. 

For the bound on a;, note that the sinh factor in the denominator of each can be estimated by 



|sinh(a + if])] = 



> 



I sinh a I 



so that 



sinh ^(a + bi^ > sinhRe ^(a + 6i)2 > |sinh(|a| - 27r)| for \a\ > 2tt. This grows exponentially 

with \a\, so it certainly dominates the polynomial growth of the numerator, and we have \ai{a + ib)\ < 1 for 
large enough \a\. □ 

Lemma 5.4. Let F{X) :— e ~^'^^'^'^^ai{X) be the integrand in h5.12\) . where x,z are fixed. If t G R™ with 
\t\ < n, then 



hi{x, z) 



F{\)d\= / F{\ + iT)d\. 



(5.20) 



Proof. Note first that F is analytic at A + i6 when |6| < tt, by the second inequality in Lemma [??^ Also, by 
Lemma 15.31 we have 

|F(A + z5)| < e-l"l'l^l/8 (5.21) 

as soon as |A| > c{9). 

We view Jj^^ F{\) dX as m iterated integrals and handle them one at a time. For 1 < k < m, suppose 
we have shown that 



F{X) dX 



F{Xi + in, . . . , Afc-i + iTfc^i, Afc, . . . , Am) dXi . . . dXr. 
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(5.22) 





na-\-iTk 


pa 




/ ^ + 


/ ^ 


/ —a 







Continuity of F and (|5.2ip show that F is integrable, so we may apply Fubini's theorem and evaluate the 
dXk integral first: 

/ F(A) dA = / . . . / F{Xi + iTi, . . . , Xk-i + iTk-i, Xk, ■ ■■ , A™) dXkdXi . . . dXm- 

Now 

/ F(Ai + in, . . . , Afc_i + iTk-i, Afc, . . . , A™) dXk 

JR 

/a 
F{Xi + iTi, . . . ,Xk-i + iTk-i,Xk, ■ ■ ■ ,Xm)dXk. 
-a 

Since Afc -F(Ai + iri, . . . , Afc_i + iTk-i, Xk, ■ ■ ■ , Xm) is analytic for |ImAfc| < (which holds because 
Kn, ■ • ■ ,Tfc)| < |t| < tt), we have 

F{...,Xk,...)dXk = 

where the contour integrals are taken along straight (horizontal or vertical) lines. But as soon as a exceeds 
c{6) from Lemma [5731 (j5.2ip gives 

\F{Xi +iTi,...,Xk-i +irfc_i,Afe,...,A„i)| dXk 

< ^^g-|x|2|(Ai,...,Afc_i -a,Afc,...,A,„)|/8 

<7re-l="l'l"l/*^0asa^oo. 
A similar argument shows the same for f", . i^, so we have 



F(Ai + iTi, . . . , Xk-i + jTfc_i, Afc, . . . , Am) dAfc 

F(Ai + in, . . . , Afc_i + iTfc_i, Afc, . . . , Am) dAfc 

-OO + lTfc 

F(Ai + iTi, . . . , Afc_i + iTfc_i, Afc + iTfc, . . . , Am) dAfc. 



Thus applying Fubini's theorem again, we have shown 



F{X) dX 



F{Xi + iTi, . . . , Afc_i + iTfc_i, Afc + iTfc, . . . , Am) dXi... dA„ 



(5.23) 
□ 



Applying this argument successively for fc = 1, 2, . . . , to establishes the lemma 

For the remainder of this section, we assume that \z\ < Bi |a;|^, so that 9 < 9o{Bi). We next show that 
the contribution from A far from the origin is negligible. 

Lemma 5.5. There exist r > and a constant C > such that 

-'l'^^+''^^^''^'^ai{X + i9z)dX < 



e 1 



B(0,r)C 



C 



(5.24) 
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Proof. From Lemma 15.31 if r > c{9o) we have 



B(0,r)C 



e 4 



dA < / e" — 1^1 dX 



B(0,r)C 



oo 



J r 



OO 



C 



where LOm-i is the hypersurface measure of 5*™ . □ 

We can now apply a steepest descent argument. As a similar argument will be used later in this paper 
(see Proposition 16. 7p . we encapsulate it in the following lemma. 

Lemma 5.6. Let S C M'"' for some k, r > 0, B{0, r) the hall of radius r in R™, and g : B{0, r) x E R, 
k : X [-r, r] x S ^ C 6e measurable. Define F : R^" x S ^ C fty 

Fix,a):^[ e-l^l'9(^''")fc(a;,A,a)dA. (5.25) 

Js(0,r) 

Suppose: 

1. There exists a positive constant bi such that ^(A, cr) > bi |A|^ for all A € B{Q,r),a € S; 

2. k is bounded, i.e. k2 :~sup,j.^g2n xeB{o,r).ae'E\^i^i'^T'^)\ < 
Then there exists a positive constant C2 such that 

\F{x,cj)\<^ (5.26) 

for all 1' > 0, (T e S. 

// additionally we have: 

3. There exists a positive constant 62 such that g{X, a) < &2 |A|^ for all A £ B{Q, r),a (£ S; 
4- There exists a function e : ~* [0, r] such that limp_>+oo p^if) = +00, and 

fci inf Rek(x,X,a) >0. (5.27) 

a;eR2i,Ae-B(0,e(|x|)),o-eE 

Then there exist positive constants C[ and xq such that for all \x\ > xq and cr € S we have 

ReFix,a)>^. (5.28) 



Proof. The upper bound is easy, since 



\F{x,a)\ <k2 e-l^l'''il^l'dA 

Js(0,r) 



k. 



-b.lM'dX 



I 1 / 

Fl J BiO,rx) 

<^ f e-^^l^l^dA 



fc2(V&l)'"/' 

\xr ■ 
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For the lower bound, let 



J B{Q,r)\B{OA\x\)) 

F2{x,a):^ [ e-l^l's(^''-)fc(x,A,a)dA 
Jb{o,4\x\)) 



so that F = Fi + F2. Now we have 



\Fi{x,a)\ < k2 / e-l"l ''il^l dX 

<k2 [ e-l-l'^^I^I'dA 

jR'"\S(0,e(|K|)) 



JR'"\B(0,|£c|€(|x|)) 



where we make the change of variables A' = |a::| A. For F2 we have 

ReF2ix,a) > h [ e-l^l'^^l^l' dA 



FI J B{Q,\x\e{\x\)) 

So we have 

|xr"ReF(x,a) > l^l™ ReF2(x, ct) - Ha;!" a)| 

>fci/ e-'^\^'Wx' ~ k2 f e-''\^'Wx' 

JB(a,\x\e{\x\)) JTSL'^\BiO,\x\e{\x\y) 

^fci(V62)"/'-0>0 
as I a; I ^ 00. So there exists so large that for all > xq, 

ReF{x,a) > lfci(V62)'"/'^ (5.29) 

as desired. □ 

We need another computation before being able to apply this lemma. 



Lemma 5.7. Re y^A+TM)^ coth y^A+TMp > 6* cot 6*, with equality ijf X = 0. 

Proof. We first note that the function /3cot/3 is strictly decreasing on [0, tt). To see this, note ^/3cot/3 = 
By Corollary [S31i^(/3) > 0. In particular, /3cot/3 < 1. 
Next we observe that for a e M, /3 G [0, tt) we have 

Re((a + iP) coth(a + il3))> f3 cot f3 (5.30) 

with equality iff a = 0. This can be seen by verifying that 

2 

cosh^ a — cos^ P 



^ ,, , , „ „ sinh a(a coth a — Scot /3) ,^ „,s 

Re((a + i/3)coth(a + i/3)) -/3cot/3 = — (5.31) 



which is a product of positive terms when a 7^ 0, since acotha > 1 > /3cot /? and cosh a > 1 > cos^ /9. 
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Therefore, we have 

Re ^/{x + iez^coth y/{\ + iezy > (im y/{\ + iezy^ cot (im ^y{x + iezy^ (5.32) 

>ecote (5.33) 



because < Im y^A + iM)^ < 9 < n by Lemma 15.2 



If equaUty holds in (|5.33p . it must be that Im \/{X + iOzY — 6. By Lemma [5.21 A and z are parallel, so 



yJ[X + iezY =±\X\+ie. if equality also holds in jSSH), we have 

Re(± I A| + i6) coth(± \X\+i6) = e cot 9 
so by (|5.30p it must be that |A| = 0. This proves the claim. □ 
Lemma 5.8. Given r > 0, there exist constants bi,b2,bg, > depending only on r and 9o such that 

hi |A|^ < ReijiX + 102,9, z) < 62 |A|^ (5.34) 

and 

|Im?/'(A + i9z, 9,z)\< 63 |A|^ (5.35) 
for all X e B(0, r) C R"\9 e [0, 9o], z e S""-i C M™. 



Proof. Note first that tp{X + i9z, 9, z) is smooth for 9 G [0, 6*0] since Ini {X + i9z) < 9 < 9o < ir, so that we 
are avoiding the singularities of wcothw. 

We have ip{i9z, 9, z) — and V \ip{i9z, 9, z) — 0. We now show the Hessian H{i9z) of at i9z is real 
and uniformly positive definite. 

By direct computation, we can find 



so that for u e 



-^{X,9,z)^A-^^')^-^ + ^'^^^is.,-^] (5.36) 



dX,dXj^^ ' ' ' ' ' A2 ^ y A2 



and in particular 



li[X)u ■ u = i^'i^tVX')^^^ + t ^ ^ ^u\' - ] (5.37) 



H{i9z)u ■ u = v'{9){z ■uf+'^ - z ■ uf^ 
= \u\' p(0) + ^(l-.) 



where s := (^j^) , so < s < 1. Note this is a real number whenever u G M™. Thus we have H{i9z)u ■ u 
written as a convex combination of two real functions of 9, so 

H(i9z)u ■ u > \u\^ min{4r-, u' (9)} > c \u\^ (5.38) 

9 

where c is the lesser of the two constants provided by Lemma [3. 121 and Corollary 13 . 1 31 respectively. This is 
valid for 9 > and hence by continuity also for 9 = 0. 
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By Taylor's theorem, this shows that (|5.34p and (|5.35p hold for small A. The upper bounds thus auto- 
matically hold for all A G B{Q,r) by continuity. To obtain the lower bound on Rei/;, it will suffice to show 
ReV' > for all A 7^ 0. But we have 



ReiP{X + iez,e,z) = eiy{0) - Re f{i0z,e,z) + Rc ^/{X + i9z)^ coth ^/{{X + iBzf 

a2 



H0) - - 

-9 cot 6* + Re 



sin^ei 



Re 



y/iX + ieFf coth y/iiX + i9z)^ 



y/{x + iez)^ coth y/iix + ieFf 



> 



by Lemma 15.71 with equality iff A = 0. 

The proof of Theorem 15.11 can now be completed. 



□ 



Proof of Theorem[5Ji We establish ([51]) first. We can apply LemmaESwith S [0, 6*0] x S"'-'^, a = {9, z), 
r the value from Lemma 15.51 and 

g{X, [9, 2)) := 7 Re V'(A + i9z, 9, z) 



k{x, A, {9, z)) e'^ lm^i^+^ez,e,z)^^(^x + i9z). 

The necessary bounds on g come from ()5.34p . For an upper bound on k, we have \k{x,X,{9, z))\ = 
|ao(A + i9z)\, which is bounded by the fact that (A, 6*, z) ranges over the bounded region 5(0, r) x [0, ^o] x 
5"'""^ which avoids the singularities of ao. 

Now for the lower bound on k. By direct computation, we have aQ{i9z) — (^j^) > 1; by continuity 

there exists S such that Ree*''ao(A + i9z) > \ for all |A| < 5 and \s\ < 6, where s e M. If |A| < ^/bs, 

where 63 is as in (|5.35p . we will have \x\^ |Im'0(A + i9z)\ < 5. Thus set e{x) := min{(5. 



-2/3 



S/b^}, so that 

Refc(x, A, (9, z) > ^ for all |A| < e(x) and all {9, z) G S, and limp^oo pe(p) = hmp^oo P^^^S/b^ = +00. 

Thus Lemma 15.61 applies, and so combining it with Lemmas 15. 41 and 15. 51 we have that there exist positive 
constants C, C{ , C2 , xq such that 



c 



:2m 



< ptix,z) < 




-2m 



(5.39) 



Then taking do = xq will 



whenever > xq. We can choose xq larger if necessary so that |a;|~ S> |x| 
establish ([5?T|) . 

For qi, the upper bound is similar; \ai\ is bounded above just like |ao|, establishing (|5.2p . 

For (|5.3p . we cannot necessarily bound both \qi\ below simultaneously, but it suffices to take them one 



at a time. For < 9{x,z) < j, we have ai{i9z) — cos9 (^j^) > so by the above logic we obtain 
the desired lower bound on \qi\ for such 9. li j < 9 < 9o, we estimate 92 in the same way, since we have 



a2{i9z) 



> 



□ 



6. Polar coordinates 

In this section, we obtain estimates forpi(a;, z) and \'Vpi{x, z)\ when \z\ > Bi \x\^, where Bi is sufficiently 
large. This means that 9{x, z) > 9q for some 9o near tt. Note that by Corollarv l3.15[ we have d{x, z) x y^jzj 
in this region. 

We first consider pi and show the following. 
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Theorem 6.1. Form odd, there exist constants Si, do such that 



or, equivalently, 



Pi(x,z) 



Pi(x,z) 



d(a;,z)2"-"-i 
1 + {\x\d{x,z)Y~^' 



(6.1) 



(6.2) 



for \z\ > Bi \x\ and \z\ > do (equivalently, d(x,z) > do). 

The effect of tlie requirement tlrat \z\ < Bi\x\^ in tlie previous section was to ensure tliat tlie critical 
point i9z stayed away from the singularities of the integrand. As Bi — > cx), the critical point approaches the 
set of singularities, and the change of contour we used is no longer effective; the constants in the estimates of 
Theorem 15 ■ 1 1 blow up. In the case of the Heisenberg groups, where the center of G has dimension m — 1, the 
singularity is a single point, and the technique used in and [H is to move the contour past the singularity 
and concentrate on the resulting residue term. For m > 1, the singularities form a large manifold and this 
technique is not easy to use directly. However, by making a change to polar coordinates, we can reduce the 
integral over R™ to one over M; this replaces the Fourier transform by the so-called Hankel transform. (A 
similar approach is used in [23i] in the context of estimates for the analytic continuation of pt.) When 
m is odd, we recover a formula very similar to that for m = 1, and the above-mentioned technique is again 
applicable. 

For the rest of this section, we assume that m is odd. 
For m > 3, we write (I4.2p in polar coordinates to obtain 



Pi{x,z) 



(27r)-™(47r)- 
(27r)-'"(47r)- 



e'P^'dae- 



dae 



- p coth p 



- p coth p 



P 



sinhp 
\ sinh p 



P 



-^dp 
-^dp 



(6.3) 
(6.4) 
a in the da 



since the integrand is an even function of p. (To see this, make the change of variables a - 
integral. It is not true when m is even.) 

The da integral can be written in terms of a Bessel function. Using spherical coordinates, we can write 
for arbitrary v G S*™"^ and w G C, 



' da — 



2n- 



sm 



' (p dtp 



(27r)"/^ 

,,,m/2-l 



(see page 79 of j21j) 



cos(u; cos ip) sin" 

Jrn/2-l{w) 

Re^H(^) ,(.) 

— l m/2— 1\ ' 



' pdip 



(by symmetry) 



where Hi,{w) is the Hankel function of the first kind, defined by H^{w) — J,y{w) + iY^{w), with Y^, the Bessel 
function of the second kind. Page 72 of 2lj has a closed-form expression for which yields 



S„,{w) - 2(2^)' 



Re 



m-1 
2 



,jm—l Z-^ 



i-iw)'' 
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(6.5) 



where the coefficients are 

(m-fc-2)! 

2— — - k)\{k - ly. 

The reason for the use of the flankel function is the appearance of the e™ factor, which gives us an 
integrand looking much hke that for pt when to = 1 . This will allow us to apply similar techniques to those 
which have been used previously for m = 1. We have 



Pi(x,z) = (Re) e'H^I-^Pcothp (gg) 

f._^ J ~oo smn p 

{m-l)/2 ^ 

E c™,fe|^l'""+'e-i'^(-^-)^ / e~'^^^P^'^^au{p)dp (6.7) 

I 1 J —oo 



k=l 

where, using similar notation as before, 

il){p, e) := -iv{e)p + p coth p 5— (6.8) 

sin 6 

(-^)"H/')' (6-9) 

The constants and coefhcients have all been absorbed into the Cm,k'i we note that ci_o > 0, c„i_fe > for 
/c > 1, and Cm.o = for to > 1. We dropped the (Re) because the imaginary part vanishes, being the 
integral of an odd function. 
For TO = 1, we can write 

Pi(x,z) = (47r)-"e-3''("'")' / e"-^'f'^P^'>^ ao{p) dp (6.10) 



The integrals appearing in the terms of the sum in ()6.7p . as well as in ()6.10p . are all susceptible to the 
same estimate, as the following theorem shows. 

Theorem 6.2. Let D G C be the strip D = {0 < Imp < 37r/2}. Suppose a{p) is a function analytic on 
D\{iTr}, with a pole of order n at p = in, a{i9) > 1 for Oq < 9 < n , and Jjj \a{p + 3i7r/2)| dp < 00. Let 



00 



h{x,z):^ e-^'f'^P''^^a{p)dp. (6.11) 

J —OQ 

There exist Bi , do such that 

Re h(x, z) X ^ p (6.12) 

l + {\x\^\r-h 

for all (x, z) with \z\ > Bi \x\^ and \z\ > dg. 



The proof of Theorem 16.21 occupies the rest of this section. Theorem 16.11 follows, since Theorem 
applies to each term of (|6.7|) (note each Ok satisfies the hypotheses), and the k = (to — l)/2 term will 
dominate for large \z\. 

An argument similar to Lemma |5.4[ using the fact that Lemma 15.31 applies for |6| < 27r, will allow us to 
move the contour to the line Imp = 37r/2, accounting for the residue at in: 



00 



h{x,z):^ / e-T^^(''+3W2,e)^(^^ 3-^/2) dp + Res(e-^'''(^'^)a(p);p = i7r). (6.13) 



hi(x,. 



The following lemma shows that hi(x,z), the integral along the horizontal line, is negligible. 
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Lemma 6.3. There exists 6q < tt and a constant C > such that for all (x, z) with 6{x, z) e [9q, tt) we 
have 

< Ce-'^t^'^''/^ (6.14) 

Proof. Observe that coth(p + 3i7r/2) = tanhp. So 

Re + 31^2,6*) = ptanhp+ — ;/(6') 



2 ^ ^ sin^fl 
Therefore we have 



< g-^(¥Ke)--fl7) I |a(p + 3j7r/2)| dp 



asrtanhr > 0. The integral in the last Une is a finite constant, since a(- + 3i7r/2) is integrable by assumption. 

However, for 9 sufficiently close to tt, we have h'{6) > ^ ^.^^ ^ . (If /?((?) i/(0) (^^jfrg^ , we have 
limei^ /?((?) = I/tt and limgi^ /3'(6') = — 2/7r^ < 0. Indeed, ^^ > 0.51 suffices.) Thus for such we have 

\hi{x, z)\ < Ce"^^ = Ce-'^("'")'/^ (6.15) 

□ 

To handle the residue term /i,., write it as 

hr{x,z)^i> e-^^('''«)/4(p)dp. (6.16) 

JdB(iiT,r) 

We can choose any r S (0, tt) because the integrand is analytic on the punctured disk. To facilitate dealing 
with the singularity aX 6 ~ -k, we adopt the parameters 



s :— TT ~ 0{x, z) 
y.^TT \xf /s. 



(6.17) 



Note that 

y/s^\z\, y^\x\y/\7\. (6.18) 

If we let (compare (|5.10l) ') 

(j){w, s) := —sip{i{TT — w), n — s) 
An 

s f (7r-s)2\ 

= — i^(7r — s){tt — w) + (tt — w) cot(7r — w) 5 (6.19) 

in \ sin s J 

F(y, s) s"-i / e-^'^('^^")a(i(7r - w)){-i) dw (6.20) 

JdB{0,r) 

we have 

hr{x,z) ^ s-^''-^^F{y,s). (6.21) 

Note we have made the change of variables p ~ i{TT — w) from (|6.16p to (|6.20p . 

Observe that F is analytic in y and s for s ^ kn, fc G Z, so we shall now consider y and s as complex 
variables. The factor of s"^^ in F was inserted to clear a pole of order n — 1 at s = 0, whose presence will 
be apparent later. 
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Computing a Laurent series for cp about (z7r,7r), which converges for < |5| < tt, < \w\ < tt, we find 

1 w s 

with U analytic for \s\ < tt, \w\ < tt. Also, by the hypotheses on a, 

a{i{Tr -w)) =w''^V{w) (6.23) 
where V is analytic for < 7r/2 and V^(0) > 0. Thus we have 

F(y,s)^s"-^ (f e-yi^-^~^-'^'^-^''^)w-''V{w)i~i)dw (6.24) 

JdB{0,r} 

The constant term in the expansion of -0 is slightly inconvenient, so let G{y, s) ~ e^/^F(y, s). Then: 
G(y, s) = s"-i / e^(3^+3?^+"^('"'"))«;-"y(w;)(-i) dw 

JdBiQ,r) 

r ^ k k 



fc! V4s Aw 

k=0 



oo t. 

2/ 



fc=0 



where we let 



(6,26) 

(^— + — + sC/(io, s) j u;""F(w)(-i)dw. (6.27) 

The interchange of sum and integral in (|6.25p is justified by Fubini's theorem, since for fixed s U{s, •) and 
V are bounded on B{0,r), and thus 



i^jB{o,r) k\ V4s Aw ^ ' ') \w 



\w\—r j y ' r 



= 27rr -+ sup |l^(u;)| exp |y| + ^ + |s| sup \\J{w,s)\\ < x>. 

We now examine more carefully the terms gu in (|6.26ffS??7|) . 
Lemma 6.4. //.gfc *s defined by {6.2T\l , then: 

1. gk is analytic for \s\ < Sq; 

2. There exists C = C{sq) > independent of k such that \gk{s)\ < for each k and all \s\ < sq; 

3. For k < n — 1, gk{s) = s"^^^^hk{s), where is analytic for \s\ < sq. In particular, gk{0) = for 
k < n — 1. 

4- For k > n — 1, 5fc(0) > when k + n is odd, and gk{0) = when k + n is even. 
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Proof. By the multinomial theorem, 

5'=(^)= E ( l (^y (^)\su{w,s)rw-vH{-i)dw 



(6.28) 



y f ^ V^^''^''^ / w'''^'"s~^''-^-"^-^sU{w,s)yV{w){-i)dw (6.30) 



a+b+c— 
a— 6— n< — 1 



since for terms with a—b — n> 0, the integrand is analytic in w and the integral vanishes. Now the integrand 
of each term of (|6.30p is clearly analytic in s, hence so is gk itself, establishing item[T] 

For itemUl let Uq := sup|„,|^^ |g|<^^ \U{w,s)\, and Vq := sup|^|^^ |l^(w)|. Then for \s\ < sq, 

\9k{s)\ < J2 (^iy-^''^'\^^r)r--'-'''so^''-''-"^-'{soUorVo 



a-\-b-\-c—k 
a— 6— n< — 1 



a+o+c— 
a—b—n<—l 

a+b+c=k \ ' ' / \ "/ 

so that a constant C can be chosen with gk (s) < , establishing item [H 
For item [HI suppose fc < n — 1 and let hk{s) — s'^~"+^5fc(s), so that 

hk{s)= y f t / w''-^-''^s-^''-^-^\sU{w,s)YV{w){-i)dw 

a—b—n<—l 

But a — b~k<a — k<0 since a < fc by definition, so only positive powers of s appear, and is analytic 
in s. 

For item 131 we see that when s = 0, each term of (|6.30p will vanish unless c — and a — b — n = —1, i.e. 
a + b ^ k and a — b — n — 1. Iffc and n have the same parity, this happens for no term, so 5fe(0) = 0. If fc 
and n have opposite parity, this forces a — {k + n ~ l)/2, b = (fc — n + l)/2, both of which are nonnegative 
integers. In this case 

9k{s)^(„ , ^ <(w-W{w){-z)dw 



{k + n- l)/2 
fc 

{k + n- l)/2 



4"''27rF(0) > 



since V(0) > 0. □ 

From this we derive corresponding properties of the function F. 
Corollary 6.5. Let F{y, s) be defined as in i6.20\) . Then for all sq < n: 

1. F is analytic for all y and all < s < sq- 
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2. We may write 



/c! 



(6.31) 



with /ifc, iJ analytic for all y and all < s < sq. Furthermore, h„-i{0) > 

3. F{y, 0) > for all y > 0. 

Proof. We prove the corresponding facts about G = e^l'^F . By items [T] and [5] of Lemma [6. 4[ we have that 
G is analytic for \s\ < sq and all y, since the sum in (|6.26p is a sum of analytic functions and converges 
uniformly. By item [3] we have that 



"~1 n-l-fe °° 



k=0 



fc=0 



{n + k) 



■gn+k{s). 



And by items El and a G(zj,0) = EZn-i 



> for aU y>0. 

Proposition 6.6. For all yi > 0, there exist S > 0, and < C( < C2 < oo such that 

G[y"~' < ReF{y,s) < \F{y,s)\ < G',y"-' 
for all < y < yi , < s < Sy. (Here we are treating y and s as real variables.) 



□ 



(6.32) 



Proof. Let if be a positive constant so large that |/ifc(s)| < K and \H{y,s)\ < K for all < ?/ < 
0<s<2/i,fc<n— 1. For any ^ < 1 and all s < dy < yi, we have 



Re G(y, s) = / Rc /i„_i(s) + V Re hkis) + Re H{y, s) 

[n — 1)1 f — i k\ 



,,1-1 



fc=0 



.n— 1 rri— 1 — fe 



'■if 



Re/tn-i(g) 
(n- 1)! 



fc=0 



^E 



fc! 



rri— 1 — fc 



y"if 



/fc! 



y^K. 



Since /i„_i(0) > 0, we may now choose 5 so small that the bracketed term is positive for all < s < 5yi. 
Then there exists yo > so small that for all < y < j/o, we have ReF(2/, s) > e~^°/^ ReG{y, s) > C[y^~^ 
for some > 0. On the other hand, 



\F{y,s)\ < \Giy,s)\ 



n-l 

^E 



k „n— 1 — 



y s 



fc! 



^E 

fc=0 



fc! 



|;ife(s)|+2/"Reii(y,s) 



Again, for small y (take yo smaller if necessary), we have \F{y,s)\ < C22/"~^. 

It remains to handle yo < y < yi- But this presents no difhculty; as F{y, 0) > for all y > 0, and F is 
continuous, there exists S so small that 

inf ReF(y,s)>0. 

yo<y<yi,o<s<5yi 



This completes the proof. 
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□ 



Proposition 6.7. There exists ?/i > 0, sq > and constants Ci, C2 > such that 



^<ReF{y,s)<\Fiy,s)\<^ (6.33) 



for all y > yi, < s < sq. 



Proof. Here the Gaussian approximation teclinique of Section [5] is again applicable. We will fix the contour 
in (|6.20p as a circle of radius r = s, parametrize it, and examine the integrand directly. Thus let w — se^^ 
in (fOIll to obtain 

g-y0(.e'-',.)^(^(^ _ se*7))se*7 (6.34) 

-TT 

We shall apply Lemma l5.6l with m = l,A = 7, r = 7r, a; = y/y. Let 

5(7,5) = Re 0(se*^,s) (6.35) 
/c(^, 7, s) = e~*v^' I'"'^('"^"^■^)s"a(^(7^ - se'^))e''^ (6.36) 

Since 0(s, s) = and = s is a critical point of (^(w, s), we have 

^0(.e-,.)|^^„ = lJ\s,s)i^sf = f!^^:^ (6.37) 

which is bounded and positive for all small s (recall v{t: — s) ^ s^^). Thus there exists SQ,e small enough 
and constants hi , 62 such that 

617' < 3(7,5) <&27' (6.38) 
for s < sq, I7I < e. Also, we have from (|6.22p that 

(/>(se'^, s) = i - i COS7 - sU{se''',s) (6.39) 

so that by taking Sq smaller if necessary, we can ensure 5(7,5) > for all s < sq Eind e < I7I < tt. Thus 
(|6.38p holds for s < sq and all 7 € [— tt, tt], with possibly different constants &i, 62- 

Boundedness of k follows from the fact that a has a pole of order n at in, so s^a{i{TT — se*^)) — V{se^'^) 

is bounded for small s. Finally, since ^^^(se*''', s)|^_p > and 1^(0) > 0, the argument used in the proof 
of Theorem 15.11 shows that the necessary lower bound on k also holds. Then an application of Lemma [5?6l 
completes the proof. □ 



Proof of Theorem ] 6. 21 Choose yi,so so that Proposition 16.71 holds, and take Bi large enough so that 
0(x,z) > TT — s when \z\ > Bi \x\^. Use this value of yi and choose a S such that Proposition 16.61 holds . and 
take do large enough that s < Sy when \z\ > do (see (|6.18l) '). So for such (x, z), either (|6.32p or (|6.33p holds; 
which one depends on the value of y — y{x, z). We can combine them to get 

n-l n-1 

" < Re^(2/,s) < |F(y,s)| < C' ^ (6.40) 

1 + 2 1 + 2/2 

Inserting this into (j6.2ip and using (|6.18p . we have (in more compact notation) 

iri-l 



hr{x,z)^(^y ^ i -. (6.41) 

^SJ l + y«-5 l + (bU/^)«-4 



By Lemma ()6.3p . hi is clearly negligible by comparison, so Theorem 16.21 is proved. □ 
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A similar argument will give us the estimates on Vpi and q2 which correspond to Theorems 14.41 and 14.51 
Theorem 6.8. Form odd, there exist constants Bi,do,C such that 



\Vpi{x,z)\ 



l + {\x\d{x,z))"+2 



and 



\q2{x,z)\<C 



1 + {\x\dix,z))"~^ 



(6.42) 



(6.43) 



whenever \z\ > Bi \x\ and d{x, z) > do. 
Proof. Applying (|4.5p to (|6.6|) . we have 
1 



Vpi{x, z) = --(27r)-'"(4^)-" \x\ (qiix, z)x + q2{x, z)J,x) 



where 



qi{x,z) = - 



2 dpi{x,z) 
R d\x\ 

(m-l)/2 



q2{x,z) 



- ^ Crn.k \z\ 

dpi{x,z) 
d\z\ 

(m-l)/2 



k—m+1 



Jp\z\-^^pcoth p 



sinh p 



n+l 



coshp)(— ip)*^ dp 



E 

A:=0 
(m-l)/2 I- 



/OO 2 / \ 

gVl.l-^pcothp (-Zp)'^-dp 
-oo Vsmhp/ 



E 

fe=0 



C-Tn^k I ^ I 



i/c— m+1 



pVkl-^^PCothp 



sinh p 



sinh p 



Each integral can be estimated by Theorem l6.2l For qi , each integral is comparable to e * d.(x,z) 



Inl- 



and the k = (m — l)/2 term dominates, so 

\qi{x,z)\ : 



|n-(m-l)/2 



i + (kK/|z|)«+t 



(6.44) 



The appearance of the extra minus sign in qi is to account for the fact that cosh(z7r) = — 1, but Theorem 
16.21 requires that a(A) be positive near A = in. 

For q2, each integral is comparable to — ^ ^ —^Q-\d(x,z) ^ ^^^^ ^^i^ k — (m — l)/2 term of the second 

sum dominates, so 



\q2{x,z)\ 



|n-l-(m-l)/2 



(6.45) 

i + {\x\^\z\r~i 

which in particular implies (|6.43p . To combine (|6.44p and (|6.45p . note that for \x\^ \z\ bounded we have 

\qi{x, z)\ X |^|"-(™-i)/2 e-|rf(-,-)^ |^2(x, z)\ X |^|"-i-(™-i)/2 ^-M-^^r (6.46) 
so that the qi term dominates, and 

(6.47) 



\Vpiix,z)\ X \x\ e-idix,zr 
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For \x\^ \z\ bounded away from we have 



|(7i(x, z)\ X |xr"~5 e-J'^^"'") 



(6.48) 



so that the qi term dominates again (""^ is bounded by assumption). Thus 



I in— (m— 1)/2 

\Vpi{x, z)\ X |a;| ^ = ^e-^'^'"'")' (6.49) 

which is equivalent to the desired estimate. □ 
7. Hadamard descent 

In this section, we obtain estimates for pi{x,z) and |Vpi(a;,z)| for \z\ > Bi \x\^ , \z\ > cLq, in the case 
where the center dimension m is even. The methods of the previous section are not directly applicable, but 
we can deduce an estimate for even m by integrating the corresponding estimate for m + 1. As discussed 
in the remark at the end of Section 21 this is valid even though there may not exist an i7-type group of 
dimension 2n + m + l with center dimension m + 1 , since the estimates we use are derived from the formula 
and hold for all values of n, m. 
We continue to assume that \z\ > Bi \x\ and \z\ > do for some sufficiently large Bi,do. To emphasize 



the dependence on the dimension, we write for the function pi in 

In order to estimate for m even, we consider We can observe that 

("■™)(x,z)= /p("'™+i)(a;,(z,z,„+i))d^™+i (7.1) 



since /jj4e*^"+i"'"+i/(A™+i)dA™+idz™+i = 2^/(0). Note that |(A,0)|r„+i = |A|a„. Now ^("'"+1) can 
be estimated by means of Theorem 16. II Using the fact that \{z, Zm+i)\ > we have that for m even, there 
exist constants Bi , do such that 

p("^"')(a;,z) X g(2«-™-2,"-^)(a;^;j) (7.2) 
whenever \z\ > Bi \x'^ and \z\ > do, where 

Jr 1 + {\x\d(X,{z,Zm+i))y 

Thus it suffices to estimate the integrated bounds given by Q'''^'^\ 
Lemma 7.1. For \z\ > Bi \x\^ and \z\ > do, we have 

^"^'"'^ i + i\x\d{x,z)r'' 



(7.4) 



We will require two preliminary computations. Since d{x,z) depends on z only through \z\, we will 
occasionally treat d as a function on M^" x [0, oo). 



Lemma 7.2. There exist positive constants ci, C2, -Bi such that for all x G M u G K with m > i?i |x| , we 



have < ci < S-d{x,u)'^ < C2 < c». 
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Proof. Let fi{e) = so that d{x,u)^ = NiVW with 6 = e{x,z) ^ y-^ 



Then 



It is easily verified that ^'(6*) > 0, v'{e) > for aU 6 G (0, tt), and ->7r>0as6'^7r. □ 
Lemma 7.3. For any a G M, there exists Ca > such that for all wq > I we have 



w^e-"" dw < Cc^w^e-""". (7.6) 
Proof. For a < 0, w" is decreasing for w > 1, so 

/"OO 

'e-'"dw<w^ / e""" = w^e-"'" (7.7) 



and this holds with Ca — 1. Now, for a nonnegative integer n, suppose the lemma holds for all a < n. Then 
ifn<a<n+l, we integrate by parts to obtain 



oo 

w^e-"" dw = w^e-""" + a I tu^-^e"'" dw < {1 + aCa^i)w^e- 



so that the lemma also holds for all a < n + 1. By induction the proof is complete. □ 



Proof of Lemma \7.1\ We make the change of variables u — |(z, Zm+i)| so that Zm+i = y — l^l^. By our 
previous abuse of notation, we can write d{x, (z, Zm+i)) — d{x,u). Thus 



i,i i + {\x\d{x,u)y y; 



V? — \z\^ 



\a+2 



y^7^y^rTRl + (k|rf(a;,zi))'3 

We used the fact that u x d{x,u)'^ where \z\ > Bi \x^ , by Corollarv l3.15l 

Now, noting that u ^ d{x,u) is an increasing function, and w i—* w^+^e^*"" is decreasing for large 
enough w, the lower bound can be obtained by 

^ ' ^""'"'-Ui ^/^\^^\i + {\x\d{x,u)r' ' ' 

- V^I^ ; ^2\z\ + ll + {\x\dix,\z\ + l))^ 

^ 1 rf(a;,|z| + l)"+^ 1^(^,1,1+1)2 

^2"^TTi + (Nrf(x,|z| + i))/5 

> „ 1 d(x,z)°+^ -i<i(a;,^)^ 

where the last line follows because u d{x,u)'^ is Lipschitz, as shown by Lemma 17.21 with a constant 
independent of x. 

Since Izj x d(x, z)^, we have that 

Q (^'^)>C. ^^^i^i^^^^^^^^e . (7.8) 
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For an upper bound, we have 



d{x, u) 



Q + 2 



Now 



< 



d{x, z) 



Q + 2 



d{x, z 



/2|I| l + 

\a + l 



e 4 



For the other term, we observe 

1 1 d(a;,u)"+2 



|,|+i v/^7~R y^TTR 1 + (Nl 



oc 



1 



a+2 



< 



y/u+ \z\ 1 + (|a:^| d{x, u))P 



< C 



|.| V2^l + {\x\d{x,u))l3 
d(a;,u)"+i 
|,| l + i\x\d{x,u))0' 



We now make the change of variables w ~ jd{x, u)^ . By the above lemma, du/dw is bounded, so 



d{x, u) 



a+l 



,|,| l + {\x\d{x,u))P 
If d{x, z) < 1/ we have 

(4u;)(a + l)/2 



idix,uf < C' 



(4u;)(" + l)/2 



1 + (2 |x| v^)/3 



where we have used Lemma 17.31 

On the other hand, when d{x, z) > 1/ \x\, we have 



^ ' ,e-'"dw <{2\x\)-^ 



+ {2\x\V^r 



id(a;,2)2 



- l + (|x|d(x,z))'3 



Combining all this, we have as desired that 



l + ilxldj^))^' 
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Corollary 7.4. Theorems \6.1\ and 1 6. (^1 also hold for m even. 

Proof. The heat kernel estimate of Theorem 16. II is immediate, given (|7.2p and Lemma l7. II 

To obtain an estimate on Vpi, we define gj"'™-* := _.A_^p("'™)(2;^ 2;), gj"'™'' '■= g^Pi"'™'' (2^1 2:), as i 
(1131). 

For qi, we simply differentiate (|7.ip to see 



m 



= / (.t, (z, z,„+i)) dz,„+i 

\2n— m+1 

" ^"^'^ _e-3^(-.-)' byLemmaO 



l + (|a;|d(a;,z))"+- 
For we again differentiate (|7.ip . Here we obtain 

(n,m) / \ f (n.m+1) f / \\ |^| , 

92 [x,z)= / '(x, (z,z„+i))-- -dz. 



|z|Q(2n-m-4,«-i) bv (lOSll 

l + (|a;|d(a;,z))"-5 



Repeating the computation from Theorem 16.81 we have the desired estimates on |Vpi| and \q2\. □ 
8. Conclusion 

An obvious extension of this result would be to obtain precise estimates for the heat kernel in more general 
nilpotent Lie groups. For step-2 nilpotent groups, a formula for the heat kernel along the lines of (|4.2p can 
be found in [5], among others. However, the additional algebraic structure enjoyed by H-type groups has 
played a major part in the analysis presented here, and its absence complicates matters considerably. A 
key difficulty is that the exponent in the formula for pt now contains expressions like Jx cot Jx , which are 
awkward to work with when Jx may not commute with its derivatives with respect to A. 
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